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ABSTRACT 


After several years of research, the image theory has been 
found acceptable to solve the sound propagation problem within 
a wedge-shaped fluid overlying either a slow or a fast bottom. 
Some further observations where done in this research. The 
primary purpose was the development and computer application 
of a new analysis, the doublet analysis, based on the image 
theory. In this analysis, the features of the image theory 
are studied from another point of view, based on a collection 
of acoustic doublets. 


iii 



TABLE OF CONTENTS 


I. INTRODUCTION... 1 

II. DEVELOPMENT. 3 

A. DEFINITIONS OF THE SYMBOLS. 3 

B. IMAGE THEORY. 6 

C. ACOUSTIC DOUBLET RADIATION WITHOUT 

REFLECTION COEFFICIENTS ... . 14 

D. ACOUSTIC DOUBLET RADIATION ENCOUNTERING 

REFLECTION COEFFICIENTS.. 17 

E. DOWNSLOPE/UPSLOPE DOUBLET ANALYSIS OF 

THE IMAGE THEORY.22 

1. General Description.22 

2. Neutral Doublet Acoustic Field.24 

3. Acoustic Field from the Pairs of 
Acoustic Doublets Ignoring the 

Bottom Reflection Coefficients....25 

4. Acoustic Field from the Pairs of 
Acoustic Doublets Encountering the 

Bottom Reflection Coefficients. . ..29 

III. COMPUTER PROGRAMS..30 

A. PROGRAM URTEXT. 30 

B. PROGRAM URTEXT 1. 32 

C. PROGRAM URTEXTIC. 32 




















SyK^NOXLiaHARV 


D. PROGRAM DOUBLTIC...32 

E. PROGRAM CALI. . . . .......34 

F. PROGRAM PLOTIFD.... 34 

IV. VERIFICATION OF THE PROBLEMS..... 35 


A. COMPARISON BETWEEN THE IMAGE THEORY 
AND JAEGER'S PE MODEL WITH RESPECT 


TO THE ACOUSTIC DOUBLET RADIATION...35 

1. Cases Tested....... 35 

2. Specific Observations. .. .35 

B. CONTRIBUTION OF THE IMAGES....36 

C. COMPARISON BETWEEN THE IMAGE THEORY 

AND THE DOUBLET ANALYSIS. . .... 37 

V. CONCLUSIONS AND RECOMMENDATIONS ....51 

APPENDIX A Tables..... 53 

APPENDIX B Program URTEXT.....58 

APPENDIX C Program DOUBLTIC. . ...67 

LIST OF REFERENCES.....76 

INITIAL DISTRIBUTION LIST. ... .77 


V 



















LIST OF TABLES 





NUMERICAL RESULTS OF COMPARISON BETWEEN THE IMAGE THEORY 


AND THE DOUBLET ANALYSIS 

1. Data for P=7°, y=4,’^, Ri=4, R2=25, Pi/Pz-^ r 

Ci=1500m/s, C2=1430m/s. 53 

2. Data for P=6°, y=2°, Ri=0.8, R2=5, Pi/p2=0.8, 

Ci=1500m/s, C2=1450m/s...54 

3. Data for p=5®, Y=4®, Ri=2, R2=400, pi/p2=0.6, 

Ci=1500m/s, C2=1410m/s. 55 

4. Data for P=4°, Y=l®, Ri=5, R2=12, Pi/p2=0.7, 

Ci=1500m/s, C2=1420m/s.56 

5. Data for p=3°, Y=1 ®a R2=10, Pi/p2=0.9, 

Ci=1500rti/s, C2=1475m/s.. .57 


VI 









LIST OF FIGOBES 


1. Image Structure for a Wedge-shaped Duct......9 

2. Illustration of the Angles ©nf 5, and 
the Angle between the n^** Image Plane and 

the Receiver.. 10 

3. The Angle of Incidence of the n*'*' Upper 

Image on the m**' Plane. ...11 

4. The Distance between Receiver and n^** 

Upper Image....12 

5. Geometry Used in Deriving the Radiation 
Characteristics of an Acoustic Doublet, 

according to Eq.23........16 

6. Geometry Used in Deriving the Radiation 

Characteristics of an Acoustic Doublet, 
according to Eq. 24..... 18 

7. Geometry Used in Deriving the Radiation 
Characteristics of an Acoustic Doublet, 

Including Reflection Coefficients, 

according to Eqs. 25 and 26..19 

8. Geometry Used in Deriving the Radiation 
Characteristics of an Acoustic Doublet, 

Including Reflection Coefficients, 

according to Eq.29.....21 

vii 











9. 


Geometry Used in Deriving the Doublet 
Analysis. .. 23 

10. Positions of the Receiver Considered in 

Program URTEXT. 31 

11. Positions of the Receiver Considered in 

Program URTEXT 1.. 33 

12. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 200m, for p=10°, 7=5.0384°, 

Pi/p2=0.9, Ci/C2=0.9. 40 

13. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 1km, for P=10°, 7=5.0384°, 

Pi/p2=0-9^ Ci/C2=0.9....41 

14. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 10km, for p=10°, 7=5.0384°, 
pi/p2=0.9, Ci/C2=0.9.. ....42 

15. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 1km, for P=5°, 7=2.5048°, 
pi/p2=0.9, Ci/C2=0.9. . ....43 

16. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 10km, for p=5°, 7=2.4952°, 

Pi/P2=0.9, Ci/C2=0.9. . . . ..44 


viii 











17. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 1km, for P=10°, Y“5.0384*’, 
pi/p2=0.94, Ci/C2»0.94 . . ......45 

18. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 10km, for P=10®, Y=5.0384*, 
pi/p2=0.94, Ci/C2=0.94..... . .46 

19. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 1km, for P=5°, Y“2.5048®, 
pi/p2=0.94, Ci/C2=0.94 . . . ......47 

20. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 10km, for Y“2.5048*, 

pi/p2=0.94, Ci/C2=0.94..... ..48 

21. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 10km, for p=10®, Y®'5*0384*, 
pi/p2=0.9999, Ci/C2=0.9999.... . 49 

22. Downslope Transmission Loss 

(dB) vs Horizontal Source-Receiver 
Distance, to 10km, for P*5*, y®2.4952®, 
pi/p2*0.9999, Ci/C2=0.9999 ..... 50 


IX 











ACKNOWLEDOIENTS 


I would like to thank Dr.A.B.Coppens for the direction and 
guidance he gave me to work on some new ideas during this 
research, and Dr.J.V.Sanders for his help in co-adyising. I 
would like also to acknowledge the support of the Computer 
Center Department and especially that of Mrs. Helen M. Davis 
and Mr. Larry Frazier. Finally, I thank my wife Evi for her 
help in preparing the write-up of this thesis. 


X 



I. INTRODUCTION 

The propagation of sound in a wedge-shaped ocean is a very 
interesting problem, which is still being investigated. A 
better understanding will be significant for scientific 
research and military applications. 

Both the image and parabolic equation (PE) models have 
been applied to this problem. Scientists who have been 
working with PE models are listed in refs. 1 and 2, and those 
with the image theory in refs. 3 through 7. At the Naval 
Postgraduate School (NPS), recent theses on this problem, are: 

1. A PE model computer program developed by Larry Ernest 
Jaeger [Ref.2]. 

2. Three other theses mostly about the image theory, which 
are listed in refs. 5 through 7. 

The most recent thesis was done by Kim Jong Rok, who, after 
a detailed development of the image theory, did a comparison 
between Jaeger's model and image theory, producing results for 
a frequency of 5Hz,a wedge angle of 10®, a source depth of 50m 
and various receiver depths. All analyses have shown that the 
image theory is acceptable. 

This research will try to further verify the reliability 
of the image theory in various cases, to identify its 
behavior, and to develop an improved analysis based on this 
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theory. More specifically this thesis will address the 
following tasks: 


1. Develop an image theory core program, which will be the 
basis of any future study. The associated computer 
program will be under the name URTEXT. 

2. Use acoustic doublet radiation to compare and analyze 
image theory results. This comparison assumes some 
matching of the water and fast bottom densities and sound 
velocities. 

3. Verify the behavior of the sound field at points where 
total internal reflections from the first images start to 
contribute. 

4. Extend a new analysis which combines the images in pairs 
of acoustic doublets, and compare with the original image 
theory. 
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II. 


DEVELOPMENT 


A. DEFINITIONS OF THE SYMBOLS 

The following symbols are used, (see Figs. 1 through 9), 
and all distances are scaled in terms of the scaling distance: 


P 

Ri 


Pi 

P2 

Cl 

Cz 

Be 

e. 

©t 

X. 

X. 


wedge angle (Eq.l) 

distance from apex to source (Eq.9) 

distance from apex to the receiver (Eq.9) 

distance along the shore between the 
receiver and the source (Eq.9) 

source angle measured downward 
from the surface (Eq.7) 

receiver angle measured downward 
from the surface (Eq.9) 

water density (Eq.l3) 

bottom density (Eq.l3) 

water sound speed (Eq.13) 

bottom sound speed (Eq.l3) 

critical angle for fast bottom (Eq.l) 

"critical angle" for slow bottom 
(defined for convenience) (Eq.4) 

transmission angle (Eq.l3) 

scaling distance for fast bottom (Eq.l) 

scaling distance for slow bottom 
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ki 

^2 

a/k2 

n 


0„ 


Run 


0 ™ 


© 


© 


Inm 


R(0un.) = 


R(0ln») = 


(defined for convenience) (Eq.4) 

wavenumber in the wedge (water) (Eq.l) 

wavenumber in the bottom (Eq.2) 

loss term in the bottom (Eq.l7) 

image number in the description of 
the image theory, and pair of acoustic 
doublets number in the description 
of the doublet analysis (Eq.7,33) 

angle of the n’'*' image measured 
from the bottom (Eq.7) 

distance between the receiver and 
the n'"*' upper image (Eq.9) 

distance between the receiver and 
the n’'*' lower image (Eq.lO) 

grazing angle of incidence on the m^*" 
plane for the n’'*' image (Eq.l3) 

grazing angle of incidence on the m*^** 
plane for the n^** upper image (Eq.ll) 

grazing angle of incidence on the m^*' 
plane for the n^*' lower image (Eq.l2) 

reflection coefficient corresponding 
to ©„^ (Eq.l8) 

reflection coefficient corresponding 
to ©!,„ (Eq,19) 

complex pressure from the upper 
family of images (Eq.l8) 
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(Eq.l9) 


complex pressure from the lower 
family of images 

distance between the images of a 
doublet 

distance of the neutral doublet 
(the source and its first surface 
reflected image) from the receiver 
receiver—neutral doublet axis angle 
defined in page 24 

distance between the receiver and 
the n*'*' upper doublet 
distance between the receiver and 
the n^*' lower doublet 
distance between the receiver and 
the n’"’’ doublet, the receiver being 
on the surface (5=0) for a given R 2 
receiver—n^*' upper doublet axis 
angle 

receiver-n^^ lower doublet axis 
angle 

receiver-n*'*' doublet axis angle, 
the receiver being on the surface 
(5=0) for a given Rj. 
as defined in eq.45 
as defined in eq.46 
as defined in eq.47 


(Eq.23) 

(Eq.31) 

(Eq.31) 

(Eq.33) 

(Eq. 34) 

(Eq.35) 

(Eq.38) 

(Eq.39) 

(Eq.40) 
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£„ = as defined in eq.48 

= as defined in eq.49 

= represent the different reflection 

coefficients of the images of an 
acoustic doublet (Eq.25) 

z s= combined reflection coefficient 


for a doxiblet (Eq.26) 

= reflection coefficients for the n^** 

upper doublet (Eq.52) 

Xn'^Vn" = reflection coefficients for the n*^** 

lower doublet (Eq.52) 

Pn'^ = complex pressure for the n^** upper 

doublet (Eq.43) 

Pn" = complex pressure for the n*"** lower 

doublet (Eq.44) 

P„ = complex pressure for the n^‘‘ pair of 

acoustic doublets (Eq.50) 


It is worth noting that the source angle y and the 
receiver angle 5 are measured downward from the surface, 
although in previous theses they were measured upward from the 
bottom. It was found more convenient to do so for the 
purposes of this research. 


B. IMAGE THEORY 

The theoretical background of the image theory has been 
presented by Coppens, Sanders, loannou and Kawamura [Ref.3]. 


6 


This section will be a brief review of the ideas and basic 
equations of the image theory. 

The assumptions are: 

1. The sound velocities and densities are constant both in 
the wedge and in the penetrable bottom. 

2. The water-air boundary is a pressure release surface. 

3. The slope is constant. 

4. All boundaries are planar. 

All distances are scaled, according to R,o«i«d = 
a fast bottom, the scaling distance is defined as the 
distance where the lowest mode attains cut off. For a slow 
bottom, the lowest mode does not attain cut off, but a scaling 
distance X, is defined for convenience in calculations. 

The definitions are 

Fast bottom 


2 sin()^ten|3 


( 1 ) 


Ir X = _ ^ _ 

^ ® 2 tanttgtanp 


( 2 ) 


cose, = ^ 


(3) 


Slow bottom 


k X - _—IE_ 

^ * 2 tan e^tanp 


(4) 
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COS0, = — (6) 

where P is the wedge angle, svibscript refers to the fluid in 
the wedge (water), and subscript 2 refers to the bottom. 

The number of images, including the source, which 
contribute to the sound field, is given by n = INT(360°/P), 
(Fig.l). The angles between each image and the bottom (Figs. 
1 through 3) are given by 

0„ = n^-y, for n odd (7) 

©„ = (n-l)P+Y, for n even (8) 

where y is the angle of the source measured from the surface. 

The distance between each image and the receiver (Figs. 2 
and 4), is given by 


= V^i"+h|-2J?,K2COs( 0„-P + 6)+ro" 
for the upper group of images, and 

= ^Rl+R^-lR^R^cos (0„+p-6) 


(9) 


( 10 ) 


for the lower group of images, where is the source-apex 
distance, Rj is the receiver-apex distance, 6 is the angle of 
the receiver from the surface, and Yq is the horizontal 
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Figure 1. Image Structure for a Wedge-shaped Duct 





Figure 2. Illustration of the Angles 0, 
n^^ Image Plane and the Receiver Plane. 





rigur« 3. The Angle of Incidence of the n^i, Upper Image on the 

Plane. 







R 


un 



2 2 2 

1 + Ra +Yo -2 Ri Rg cos (8„- p + 5) 


Figure 4. The Distance R^n between Receiver and n’^*' Upper Image. 
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source-receiver distance along the shore. All distances are 
scaled. 

The incident angle of the n’"'* image to the m^^* plane (Figs. 
1 and 3), is given by 


. , KjSin (©„-2mP)+jR sin[ (2m-l) P+6] (11) 

= —--^-_,m=l,2,... 


for the upper group of images, and 


. , H^sin (©„-2jnp)+jR 2 Sin[ (2jn+l) p-6] (12) 

sin(©i^) = -J:-1!--, m=0,l,... 

•‘'in 


for the lower group of images. 

The reflection coefficients are assumed to be the plane 
wave Rayleigh reflection coefficients, a good approximation 
for the problem. The equations for them are: 


Without atbsorption 


-^sin(©^)-^sin(©,) 

R(©^) = Pi--- (13) 

^sin(©^)+^sin(©^) 

Pi ^2 



where ©^ is the incident grazing angle of the n^** image on the 
m'^*' plane (Fig. 3) and ©t is the corresponding transmission 
grazing angle. 
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jb = 2(c,/c2)2a/A:2 (17) 

where cc is the absorption coefficient. 

The total sound field pressure contribution from the upper 
images at each point in the wedge, is given by 

= £-J^exp(-j*:iR„„) (-l)<"/2)nR(e„^) (18) 

nR(©u™) * 1, for n = 1,2 and M = (n-l)/2, 
and from the lower images, 

= E (-1) (19) 

The total pressure will be, P = P„+Pi. 

C. ACOUSTIC DOUBLET RADIATION WITHOUT REFLECTION COEFFICIENTS 
An acoustic doublet consists of two simple sources of 
equal strengths, separated by a distance d, vibrating with the 
same freqpiency but 180® out of phase with each other [Ref.8] . 

Assiime as a reference that the source closer to the 
receiver is the negative one and the other one is the positive 
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(Fig.5). The pressure at point (r,6) due to the positive 
source is 


A g il®t-A:(r+Ar,) ] 

r+Ar, 


and that due to the negative source is 





( 21 ) 


where r is the distance from the field point to the midpoint 
between the sources, Ar^ and Arj are as shown, and the opposite 
signs account for the phase difference of the sources. The 
acoustic doublet pressure at (r,0) is then 



The process of replacing Ar^ and Arj by r and 0 is 
straightforward but tedious, and the resulting expression for 
p(r,0,t) is too complicated for efficient analysis. However, 
in the most frequently encountered cases, as it happens in the 
wedge-shaped ocean problem, the observation of the pressure is 
made at distances great compared to the separation of the 
sources. Therefore, it will be useful to derive the form of 
eq.22 suitable in the limit r»d. This is referred to as the 
far-field approximation. If r»d, Ari=Ar 2 « (d/2) sin0 and 
Ari/r«l, so that eq.22 becomes 

p(r,0, t) = -j^sin (.^Jcdsin0)e^<®‘'*'> (23) 
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Figure 5. Geometry Used in Deriving the Radiation Characteristics of an 
Acoustic Doublet, according to Eq.23. 
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If the signs of the sources are reversed (Fig. 6), then 
eq.23 changes sign too, and becomes 

p(r, e,t) = j^sin (.i.kdsine)e^<®‘'*'^> (24) 

r 2 


D. ACOUSTIC DOUBLET RADIATION ENCOUNTERING REFLECTION 
COEFFICIENTS 


Consider again the situation of Fig.5, which was analyzed 
in the previous section. Suppose now that the free acoustic 
fields of the positive and negative sources encounter 
reflection coefficients % and V respectively. Those 
reflection coefficients are assigned to correspond to grazing 
angles 0^ as shown in Fig. 7. Then eq.22 becomes 



(X 






e _lir ® 

iTKFJr ^l-A£Jr^ 


(25) 


Combination of eq.23 and eq.25 gives 


A 

~r 


-j^ Maine 


-V 


j^kdmLnd 


-j2^z sin (.^/cdsin0) (26) 


for the far field approximation (r»d) , where 
Ari«Ar 2 = (1/2)dsin0 and (Ari/r)«l, (Ar 2 /r)«l. The coefficient 
z. is the appropriate reflection coefficient for the acoustic 
doublet field. After some mathematical manipulation of eq.26, 


z = ^[(X+V)-j(V-X)cot(^*dsin0)] (27) 
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Figure 6. Geometry Used in Deriving the Radiation Characteristics of an 
Acoustic Doublet according to Eq.24. 









7. Geometry Used in Deriving the Radiation Characteristics of an 
Acoustic Doublet, Including Reflection Coefficients, according to Bqs. 
25 and 26. 
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Suppose now that the polarity of the sources changes 
(Fig. 8), and consider eq. 24. In that case, the signs in both 
sides of eq.26 change, and the solution of eq.27 for the 
doublet reflection coefficient ^ remains the same. 

Surprisingly, this reflection coefficient didn't turn out 
to be too complicated to apply in the doublet analysis of the 
image theory, because it is expressed as a function of the sum 
and the difference of the reflection coefficients of the 
individual images in the doublet. 

In this research, the choice for the difference was the 
reflection coefficient of the image closer to the receiver 
minus the reflection coefficient of the image further from the 
receiver. That can be seen from eq.27 in connection with Figs 
7 and 8. 

Combination of eq.27 with eqs. 23, 24 and 26, gives for 
the two different polarities the total acoustic doublet fields 

p' = --i [ (V“X) (.iicdsinG) +j (X+V) sin (.^.fcdsinO) ] 

p* {\\f-X) cos i-^kdiSinQ) +j ix+y\f) sini^kdsinO) (29) 

with A has been set equal to one. 

These two equations were successfully applied in the 
doublet analysis of the image theory as will be discussed 
later. 
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Figure 8 . Geometry Used in Deriving the Radiation Characteristics of an 
Acoustic Doublet, Including Reflection Coefficients, according to Eq.29. 






E. DOfVNSLOPE/UPSLOPE DOUBLET ANALYSIS OF THE IMAGE THEORY 

1. General Description 

The objective of this analysis is to combine the 
images in pairs of acoustic doublets and thus simplify the 
image theory [Ref.9]. 

The procedure starts (Fig.9) by first considering the 
1’^ and 2"** upper images, the source and its first surface 
reflected image, as a neutral doublet. It is obvious that 
this neutral doublet (n=0) does not encounter any reflection 
coefficients. 

The 3’^'^ and 4*^*' upper images (Fig. 9) constitute the 1*’' 
upper doublet, and the 1**^ and 2"“^ lower images constitute the 
I**" lower doublet. Both doublets constitute the I"’" pair of 
acoustic doublets (n=l). This procedure continues until all 
the images are organized in pairs of acoustic doublets. 

The total number of pairs of acoustic doublets is 



It is easily seen that in the worst case three images 
might be missed. In that case, those will be images at the 
end opposite to the source. Their contribution will be 
negligible, because many reflection coefficients will 
contribute to them, and among those reflection coefficients 
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will be some with very small values. Hence, those images will 
be ignored. 

2. Neutral Dotablet Acoustic Field 

The neutral doiiblet, which consists of the source and 
its first surface reflected image (Fig.9), does not encounter 
any reflection coefficients. 

Application of eq.24 gives the pressure from the 
neutral doublet, 

Pq = + j_£-sin (-iArj^dsincTo) (31) 

■^0 ^ 

sinGo = —sxnd = —b 

.iA:j^dsin0o= 

2 


If we define 


ro= 


•fci 




then the pressure has the simple form 

Po= + j-1 sin (FoYS) 

■^0 


(32) 


where Pq is the neutral doublet complex pressure, ro is its 
distance from the receiver, and 0o is the receiver-doublet- 
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:^sin(2np) 

^nO 

= 5 ^ [cos (2i2p) -^sin^ (2np) ] 

■*^no r„o 

= r„o+.^i^6 [sin (2nP) +i8cos (2np) (2np) ] 

•^«o 2 2 j.2^ 

r„' = r„o+.^i.--.S [-sin (2np) +.iScos (2nP) --jLs ^^^sin^ (2np) 

^no 2 2 

Receiver-doublet-doublet axis angles are 

sin(a"^) = _^sin(2nP+8) (38) 

for the upper doublet and 

sin(a„') = —isin (2nP-6) (39) 

for the lower doublet. 

When the receiver is at the surface this angle becomes 
the same for both doublets, 

sin(a„o) = -^sin(2np) (40) 

^nO 

Eqs. 38 and 39 to the I**' order approximation of the 
Taylor's series expansion around 6=0 are 


^ d^r„- 

"W 8,0 6=0 


dr„- 

5=0 
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dr. 


= r„6sin(2np) = 


5«0 


( 47 ) 


dsinc. 


= r„YScos(2np) = i^k^d) 6—^ 


8=0 


(48) 


? = ^ l;,a2 

-yy'^-TT- -vr-rr-f - 


d^r 


8=0 


(49) 


The above approximations and the definitions in eqs. 
45 through 49, give the pressure field for the n^** pair of 
acoustic doublets (the sum of the pressure fields of the 
corresponding upper and lower doublets) 


Pn ^ = (-I)'’!—[e'^’-sin (i>„+e„)-e^'^-sin (i^„-e„) ] 

(50) 

More specifically, the following approximations where 
applied in eq.50: 

1. r„o used in amplitude {2/r^Q)f instead of r^"^ and r„". 

2. r„ correct to the 2"*^ order approximation (4„) in phase. 

3. (Tn correct to the order approximation (£„) in 

amplitude. 

The total field is the sum over all pairs of acoustic 
doublets, plus the field of the neutral doublet: 

P = [P/fP„-] (51) 

n-0 
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4. Acoustic Field from the Pairs of Acoustic Doublets 

Encountering the Bottom Reflection Coefficients. 

The reflection coefficient of eq.27 is applied to 
eq.50. Hence, according to eqs. 28 and 29, the pressure field 
including the bottom reflection coefficients is 

i (-1)" <cos(i^„+e„)) (Y/-X/)(sin (i&„+e„)) 

^nO 

(52) 

(cos (i»„-e„)) +J(sin (i»„-ej ) (X„ +V„-) ]) 

where Vn” reflection coefficients of the upper and 

lower doublet images closer to the receiver, and , Xn 
those further from the receiver. The corresponding grazing 
angles 6^", 0y^, 0y" are illustrated in Fig. 9. 

The total field is given by eq.51, where P„=P„''^+P„" is 
given by eq.52. 
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III . COMPUTER PROGRJ^S 

The following computer programs, were created and used 
during this research. They were written in the FORTRAN 77 
computer language, on the IBM-3033 main frame NPS computer 
system. Some of the programs use plotting routines, using the 
DISSPLA graphics package, while URTEXT and DOUBLTIC do not. 
Copies of all programs are available in professors Coppens 
main frame account. 

A. PROGRAM URTEXT 

In the beginning, the program URTEXT was created to 
provide a well documented, immutable program of the wedge- 
shaped ocean image theory. The follow-on programs are 
modifications of URTEXT. 

In URTEXT, the receiver, at distance Rj from the apex, 
moves from the surface to the bottom of the wedge, on a 
circular arc (Fig.10). This program provides two distinct 
results of the image theory: 

1. Pressure amplitude vs receiver angle. 

2. Phase vs receiver angle. 

The receiver angle is measured downward from the surface. 
This program accepts input data and also gives the output data 
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Figure 10. Positions of the Receiver Considered in Program 


^APEX 


URTEXT 


3 ] 



with respect to the scaling distance, independent of the 
frequency. Appendix "B" is a copy of URTEXT. 

B. PROGRAM URTEXT1 

In this program the receiver moves horizontally. Hence, 
its depth from the surface is kept constant (Fig.11). This 
program provides two distinct results for comparison: 

1. Transmission loss (TL) vs source-receiver horizontal 
range, using the complete image theory. 

2. Transmission loss (TL) vs source-receiver horizontal 
range, using only the neutral doublet acoustic radiation. 

The inputs and outputs are frequency dependent and 
distances are expressed in meters or kilometers. 

C. PROGRAM URTEXTIC 

In this program, the receiver moves as does in URTEXT and 
it provides the same results. The only difference is the 
addition of a plotting routine, to present the data. 

D. PROGRAM DOOBLTIC 

In this program, the receiver moves the same as in URTEXT, 
and it also provides the same results. The calculations are 
done according to the doublet analysis presented previously, 
without including any absorption in the bottom (a/k 2 = 0 ) or in 
the water. The reflection coefficient used, was that of 
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equation 13. It works for a slow bottom only. It was used 
with their real reflection coefficients, this program was used 
to simplify the initial verification of the doublet analysis. 
Appendix "C" is a copy of this program. 

E. PROGRAM CALI 

This program calculates the cosine of the incident grazing 
angle of the n’^*’ image on the m^** plane, as a function of 
range. It was used to determine at what range from the source 
the total reflection contribution from the first images 
starts, as a function of the critical angle, defined as 

COS0o=Ci/C2. 

F. PROGRAM PLOTIFD 

This is a plotting program developed by Jaeger to plot his 
PE model/s data. It was modified to overplot URTEXTl data for 
comparison. 
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IV. VERIFICATION OF THE PROBLEMS 

A. COMPARISON BETWEEN THE IMAGE THEORY AND JAEGER^S PE MODEL 
WITH RESPECT TO THE ACOUSTIC DOUBLET RADIATION. 

This was done by using the URTEXTl, JAEGER and PLOTIFD 
programs. 

1. Cases Tested 

The cases tested (Figs. 12 to 22) were: 

1. Wedge angles p = 10° and p = 5°. 

2. Frequency f = 50 Hz. 

3. Source depth = 50 m. 

4. Source angles y = 2.4952°, 2.5048°, 5.0384°. 


5. 

Receiver depth = 

75 m. 


6 . 

D1 = pi/p2 = 0.9, 

0.94, 

0.9999. 

7 . 

CC = Ci/Cj = 0.9, 

0.94, 

0.9999. 


8. Source-receiver horizontal ranges from 0-1 km and from 0- 
10 km, downslope. 

In these figures, the dashed line represents the image 
theory data, the dotted line represents the dotablet (source 
and first surface reflected image) data, and the solid line 
represents the Jaeger's PE model data. 

2. Specific Observations 

1. In close distances from 0-200 m (Figs. 12 and 13), the 
image theory is almost in agreement with the acoustic 
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2 . 


doiiblet radiation curve, while the PE model does not 
agree very well. 

At larger distances, it is hard to say if the image 
theory or PE model is better, because the contribution of 
the images becomes significant. More specifically: 

a. The two TL curves seem to agree with each other, 
being far away from that of the acoustic doublet. 

b. In almost identical matching, Di=0.9999, CC=0.9999 

and at P=10° (Fig.21), all three curves were 
perfectly coincident. The same case at {5=5° 
(Fig.22), gave a difference between the image theory 
and the acoustic doublet, starting from about 3 km, 
and reaching about 2.5 db at 10.0 ]ati, but gave a good 
coincidence up to 3 km. This might have happened due 
to the more significant contribution of more images 
at larger distances. 

3. For P=10° and 0^=0.9, CC=0.9 or closer to unity, for 

short ranges up to 100 m, the effects of the bottom can 
be ignored with very small error. 

B. CONTRIBUTION OF THE IMAGES 

By using the program CALI, some calculations were made of 
the grazing angle of the n^** image with the m^** plane versus 
range. When this angle is equal to the 0^, this is the range 
where total internal reflection of the image starts. This 
test was made for the 1"'" and 2"'* lower images and the 3^“^ upper 
image, because these contribute most significantly. The 
following observations, were done for the cases considered 
previously: 


1. At P=10°, CC=0.9, the I”*" lower image contributes total 
reflection at 236 m, the 2"** lower image at 533 m, and the 
3'"'* upper image at 764 m. Almost at those locations three 
characteristic bumps appear on the image theory TL curve 
for that case (Fig.l3). On the PE model curve the first 
bump appears correctly, the second did not appear, and 
the third appeares with a difference of 33 m (Fig.13). 
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2. The same case was tested at P=5°, CC=0.9, and the 

distances obtained where 187 m,427 m, 581 m. The image 
theory TL curve'^s bumps correctly appeared almost at 
those distances (Fig.15). 


All tests done above show that the image theory behaves as 
expected. Hence, it is reconfirmed that it is deemed worthy 
to be applied in the wedge-shaped ocean problem. 


C. COMPARISON BETWEEN THE IMAGE THEORY AND THE DOUBLET 
ANALYSIS 

This comparison was done using the results of URTEXTIC and 
DOUBLTIC. Some of the results obtained can be seen in Tables 
1 to 5, in Appendix "A". Those tables clearly show the 
excellent agreement between the two programs. Many other 
cases were tested to verify those results, but are not 
included in this thesis. All results that were tested were in 
agreement. More specifically: 

1. In all cases tested, there was an excellent agreement in 
amplitude, as it can be seen from some data collected in 
Tables 1 to 5. 

2. Although very good agreement (3’^'^ decimal) is still valid 
in amplitude when the source gets close to the bottom, 
in that case appears a phase difference not larger than 
30° in most cases tested. An example can be seen in the 
data collected in Table 1. The reduced accuracy of the 
doublet approximation (Eqs. 23 and 24) , due to the larger 
d, is not responsible, because the phase difference does 
not disappear even when R 2 -Ri»d. It was seen that with 
higher order approximations of the Taylor's series 
expansion of r„, a„ in both phase and amplitude of eq.52, 
that phase difference reduces. Probably, it is not worth 
complicating the analysis by including these higher 
orders, because the accuracy as it exists will be 
acceptable in most practical cases. Also note that only 
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the amplitude is necessary for the calculation of the 
transmission loss, and that is obtained with excellent 
accuracy. 

Because of the good agreement described previously, the 
doublet analysis can be used as a substitute of the full image 
theory. 

During this research, efforts were made to simplify the 
doublet analysis. It was found that the approximations 
described in sect. E of the development, can not be simplified 
to give correct results everywhere in the wedge. 

Consider the grazing angles formed between r„* and bottom 
planes for the upper doviblets, and between r„~ and bottom 
planes for the lower doublets (Fig.9). A further 
simplification of the doioblet reflection coefficients could be 
accomplished by considering that the reflection coefficients 
correspond to these grazing angles. This method gives correct 
results only under the follow conditions; 

1. Large wedge angles (P>5®), which usually are not met in 
practice. 

2. The receiver shallower than the surface-bottom middle 
depth. 

3. Ri/R 2«1 cannot hold for the downslope case, nor Ri/R 2»1, 
for the upslope case. 

It was not possible to determine an exact threshold to 
apply such an approximation, because the total field is 
calculated encountering multiple factors. For that reason, 
the reflection coefficient as given in eq.27 was used. Some 
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explanations of the necessity to use eq.27 for the reflection 
coefficient, can be derived by examining eq.28. This equation 
is reproduced here without the phase factor ; 

p~ = -—[ (\(r-x) cos (—icdsinS)+j (x) sin (-|icdsin0) ] 

X" 2 2 

X+y is always greater than ^or real reflection 
coefficients. The smaller the difference between the 
reflection coefficients Xf¥ images in the doublet, the 
more dominant the imaginary part of the above equation will 
be, with respect to X^V* 

Consider now, small wedge angles (p<5°) , along with Ri«R 2 
for downslope and Ri»R 2 for upslope. Due to the small wedge 
angle, the lower images contribute significantly. Then, for 
these lower images, cos [ (1/2) kdsinO]—and sin[ (1/2) kdsinG]—)0, 
because of the small value of the receiver-doublet axis angle, 
under the above conditions. Therefore, in that case, the real 
part will be more dominant. 

It becomes clear from the above discussion, that in many 
cases, neither the real or the imaginary parts of eq.28 can be 
ignored. Thus, eq.28 is the correct equation to apply for the 
wedge. 


39 





rigur* 12. Downslope Transmission Loss (dB) vs Horizontal Source- 
Receiver Distance, to 200m, for p-10*, 7-5.0384*, pj/p,-0.9, Cj/C,-0.9. 
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rigurs 14. Downslope Transmission Loss <dB) vs Horizontal Source- 


Receiver Distance, to 10km, for p-10®, y-5.0384®, pi/p2-0.9, Ci/Cj-0.9. 
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Flgur* 16. Downslope Transmission Loss (dB) vs Horizontal Source- 
Receiver Distance, to 103cm, for P-5*, Y-2.4952*, Pi/pj-0.9, Ci/C2“0.9. 











17. Downslope Transmission Loss (dB) vs Horizontal Source- 
Receiver Distance, to 1km, for P*10*, Y*5.0384®, Pi/p2"0.94, Cx/C2*0.94. 
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Flgur* 18. Dovnslope Transmission Loss (dB) vs Horizontal Source- 
Receiver Distance, to 10km, for p-10*, Y“5.0384*, pi/p 2 - 0 . 94 , C 1 /C 2 -O. 94 . 
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rigur« 20. Downslope Transmission Loss (dB) vs Horizontal Source- 


Receiver Distance^ to 10km, for p-5®, y-2.5048^, pi/p 2 « 0 . 94 , C 1 /C 2 -O. 94 . 
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rigor* 21 . Downslope Transmission Loss (dB) vs Horizontal Source- 
Receiver Distance, to lOkm, for P«10*, y* 5.0384*, pt/p,“0.9999, 
Ci/C,-0.9999. 
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rigur* 22. Downslope Transmission Loss (dB) vs Horizontal Source- 
Receiver Distance, to lOlcm, for P»*5®, Y-2.4952®, p^/pj-0.9999, 
Ci/Cj-0.9999, 
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V. CONCLUSIONS AND RECOMMENDATIONS 

Image theory is commonly accepted as a solution to the 
wedge-shaped ocean acoustic propagation problem, under the 
assumptions defined in sect. B of the development. Doublet 
analysis is a successful substitute for the image theory. 
This is important for two reasons: 


1. It is a more compact solution. 

2. The corresponding computer program execution time is 
faster than that of the image theory, because it requires 
fewer loops to calculate. Some tests which were done 
have shown that t^^O. 55ti, where t^ is the doublet 
analysis computer program execution time and tj is the 
image theory analysis computer program execution time. 


It is recommended that the doublet analysis be extended 
for the cross slope case, and that the corresponding computer 
program be extended for a fast bottom and to include 
absorption in the bottom. 

It also appears possible an additional simplification, can 
be made according to the dipole approximation [Ref.8], where 
the doublet acoustic pressure without including any reflection 
coefficients 

p(r,0,fc) = -j.^sin (-iiCj^dsinS) 
further reduces to 
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p(r, 0, t) = 

r 

when kid=2k;iRiY«l. Using equations 1 and 4, this can also be 
expressed with respect to the scaling distance, as 



< 1 


for a fast bottom, and 



for a slow bottom. It is obvious from the last two equations 
that the dipole approximation can be applied whenever Ri«X. 
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APPENDIX "A" 

NUMERICAL RESULTS OF COMPARISON BETWEEN THE IMAGE THEORY AND 
THE DOUBLET ANALYSIS 

TABLE 1 


Data for 6 = 7 ° , y= 4<> .R,=4,R,=25, P,/o..=0.7, c-.=1500m/s, c,=1430m/3 


RECEIVER 

ANGLE 

PRESSURE 

AMPLITUDE 

PHASE 

ANGLE 

IMAGE 

THEORY 

DOUBLET 

ANALYSIS 

IMAGE 

THEORY 

DOUBLET 

ANALYSIS 

o 

o 

0 

0.005314 

0.005282 

58.0° 

85.4° 

o 

o 

0 

0.055618 

0.055174 

56.9° 

84.4° 

1.40° 

0.124116 

0.123650 

54.4° 

81.8° 

2 . 10 ° 

0.212473 

0.212157 

51.3° 

78.6° 

2.80° 

0.311260 

0.311706 

48.7° 

75.8° 

3.50° 

0.394972 

0.396090 

i 

46.8° 

73.7° 

O 

0 

0.434343 

0.435430 

45.4° 

72.3° 

4.90° 

0.412928 

! 

0.413046 

44.1° 

70.9° 

5.60° 

0.331062 

0.329867 

41.6° 

68.3° 

6.30° 

0.200917 

0.198261 

33.5° 

59.6° 

7.00° 

0.069522 

0.072188 

-25.4° 

-4.0° 
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TABLE 2 


Data for 3=6° , Y=2° , R,=0.8, R,=5, On/O,=0.8. c,=1500m/s. c-,=145Qm/fl 


RECEIVER 

PRESSURE 

AMPLITUDE 

ANGLE 

IMAGE 

THEORY 

DOUBLET 

ANALYSIS 

O 

O 

o 

0.000372 

0.000371 

0.60° 

0.003728 

0.003667 

1.20° 

0.007214 

0.007091 

1.80° 

0.010212 

0.010050 

0 

o 

CM 

0.012559 

0.012363 

3.00° 

0.014134 

0.013921 

3.60° 

0.014946 

0.014720 

4.20° 

0.015168 

0.014929 

00 

o 

0 

0.015224 

0.014966 

5.40° 

0.015831 

0.015523 

6.00° 

0.017739 

0.017334 


PHASE ANGLE 


IMAGE DOUBLET 


THEORY ANALYSIS 


-50.2° -46.0° 


-48.9° 


-55.3° 


-60.2° 


-67.1° 


-77.0' 


-90.5° 


-46.4° 


-49.9° -47.6° 


-52.1' 




-53.0° 


-57.8° 


-64.7° 


-74.5° 


- 88 . 1 ° 


-105.4° 


-123.9° 
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TABLE 3 


Data for B=5°. Y=4°,R,-2,R,-400 


0.6. c,g l500m/3, c,=1410m/s 


RECEIVER 

PRESSURE 

AMPLITUDE 

ANGLE 

IMAGE 

THEORY 

DOUBLET 

ANALYSIS 

o 

• 

o 

(J1 

o 

0.000087 

0.000083 

0.50® 

0.000818 

0.000817 

1 .00® 

0.001532 

0.001558 

1.50® 

0.002151 

0.002158 

to 

O 

o 

0 

0.002537 

■ ■ 1 

0.002559 

2.50® 

— 

0.002704 

0.002725 

3.00® 

0.002639 

0.002635 

3.50® 

0.002276 

0.002292 

0 

§ 

0.001699 

0.001715 

o 

O 

in 

0.000926 

0.000942 

Cn 

o 

o 

0 

0.000025 

0.000032 


PHASE ANGLE 


IMAGE 


56.3< 


62.6^ 


62.2'> 


63.0® 


62.6® 


63.3® 


63.0® 


62.4® 


63.0® 


60.5® 


-30.6® 


DOUBLET 


THEORY ANALYSIS 


78.8® 


78.8® 


78.9® 


79.0® 


wm 


79.4® 


79.5® 


_ 


78.8® 


39.8® 
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TABLE 4 


Data for B=4 


Ri=5,R,=12, Di/D,=0.7, Ci=1500ni/s, c,=1420m/s 


RECEIVER 

PRESSURE 

AMPLITUDE 

PHASE 

ANGLE 

ANGLE 

IMAGE 

THEORY 

DOUBLET 

ANALYSIS 

IMAGE 

THEORY 

DOUBLET 

ANALYSIS 

o 

o 

0 

0.040927 

0.040984 

-149.3° 

-144.2° 

O 

O 

0 

0.394236 

0.394086 

-148.6° 

-143.5° 

O 

CD 

O 

o 

0.701851 

0.701115 

-146.5° 

-141.4° 

1.20° 

0.873027 

0.871301 

-142.3° 

-137.3° 

1.60° 

0.915180 

0.912006 

-135.2° 

-130.2° 

0 

O 

o 

CM 

0.887048 

0.882283 

-124.6° 

-119.6° 

0 

O 

CM 

0.851987 

0.846756 

-112.8° 

-107.7° 

to 

00 

O 

0 

0.810276 

0.805379 

-104.5° 

-99.4° 

0 

O 

CM 

CO 

0.704744 

0.699685 

-103.6° 

-98.7° 

3.60° 

0.503485 

0.498563 

-115.2° 

-110.9° 

o 

o 

0 

0.315934 

0.320342 

-161.8° 

-158.9° 



56 
































































TABLE 5 


Data for 6=3* 


R,=1,R,=10, On/O.>=0 . 9, c,=1500m/s, c.»=1475m/s 


RECEIVER 

PRESSURE 

AMPLITUDE 

ANGLE 

IMAGE 

THEORY 

DOUBLET 

ANALYSIS 

0.03® 

0.000123 

0.000127 

0.30® 

0.001255 

0.001257 

o 

o 

o 

0.002415 

0.002414 

0.90® 

0.003352 

0.003377 

1 .20® 

0.004028 

0.004066 

1.50® 

0.004384 

0.004421 

1.80® 

0.004369 

0.004409 

2 .10® 

0.004010 

0.004030 

2.40® 

0.003339 ; 

0.003329 

2.70® 

0.002471 ! 

0.002444 

3.00® 

0.001898 

1 

0.001799 


PHASE ANGLE 


IMAGE DOUBLET 


THEORY ANALYSIS 


-54.4*’ 


-49.8*> 


-48.8* 


-49.4*’ 


-48.5® 


-47.5® 


-45.7® 


-42.1® 


-35.2® 


- 20 . 8 ® 


12.7® 



I 


-47.2® 


-47.0® 


-46.6 


■ 




-45.1' 


-43.3® 


-40.2® 


-34.1® 


- 20 . 1 ® 
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APPENDIX "B" 

******* Dickie** i(i< 

PROGRAM URTEXT 

IMAGE THEORY—CROSS SLOPE PRESSURE VS RECEIVER ANGLE(D) 
*************************************************************** 

INTEGER I,I1,N,S1,S2,N1,A 

REAL*4 B,CC,C2,D1,D2,G,PI,P1,P2,Q1,R1,R2,T,T4, 

* T6,Y0,Y1,Y2,Z1,Z2;Z3,Z4,Z5,Z6,T1 (900), 

* R8 (900) ,R9(900) ,S(900) ,C(900) ,E(900) , 

* F(900),Z,R3,AL,PZ,D 

DOUBLE PRECISION Y,W0,W1 

REAL TQQ,TQQ1,TQQ2,TQQ3 

PI =ACOS(-1.0D00) 

WRITE(*,*)'INPUT B,G,D1,CC,R1,R2,AL,Y0,A' 

READ (6, *)B,G,D1, CC,R1,R2,AL, Y0,A 
WRITE (6,2) B,G,D1,CC,R1,R2,AL, Y0,A 
2 FORMAT('B = ',F4.1,/,'G = ',F4.1,/,'D1 = ^ 

* ,F4.1,/,/CC = ',F4.1,/,'R1 = ',F4.1,/,'R2 =',F4.1,/,'AL 

* = ' ,F6.3, */,'YO = ',F4.1,/,'A = '’,14) 

i(i< i< if if Ik Tk ^ if i(i< •k it -k-k 9( i< fcii * ic i(i< if 

C INPUT PARAMETERS 

C****************************** ***************************** 
C B = WEDGE ANGLE (DEG) 

C G = SOURCE ANGLE (DEG) 
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C D = RECEIVER ANGLE (DEG) 

C Nl= # OF IMAGE POINTS = INT(180/B) 

C R1 = SOURCE DISTANCE (IN DUMP DISTANCES) FROM APEX 

C R2 = RECEIVER DISTANCE (IN DUMP DISTANCES) FROM APEX 

C YO = DISTANCE (IN DUMP DISTANCES) ALONG APEX 

C D1 = RH01/RH02 

C CC = C1/C2 

C AL = ALPHA/K2 

C A = # OF RECEIVER POSITIONS 

Q ----^-----— 

C CHOOSE SLOW OR FAST BOTTOM BY VALUE OF SPEED RATIO 
C CC = C1/C2 . 

Q -----^---------- 

**Vt**i>cyc>ic********:«:*****:*c****:*:**************************^^ 

C MAIN PROGRAM 

Qyc*^^ycV^Vtyc*****:it**i'c*yc*i*f****^"cyc**i*c:«c*:*f**?^****:»cT*f*i«f*i«r:^*****i*f******* 

N1 = INT(180./B) 

T6 = 180./PI 

B = B/T6 

G = B-(G/T6) 

C2 = CC**2 

TQQ = TAN(B) 

C__--^---—------ 

C DECISION ABOUT SLOW OR FAST BOTTOM 
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C T4 = PI/(2*SIN(AC0S(CC))*TAN(B)) FOR FAST BOTTOM 

C T4 = PI/(2*TAN(ACOS(1/CC))*TAN(B)) FOR SLOW BOTTOM 

C—-^---______------ 

c 

IF (CC.LT.l) THEN 
TQQl =ACOS(CC) 

TQQ2 = SIN(TQQl) 

ELSE 

TQQl = ACOS(l/CC) 

TQQ2 = TAN(TQQl) 

ENDIF 

TQQ3 = 2.*TQQ2*TQQ 

T4 = PI/TQQ3 

Q1 = 1/DSQRT(2.0D00) 

C 

DO 20 Ml = 1, A 
D = B-(B/A)*M1 
D2 = (Y0*Y0) + (R1*R1) + (R2*R2) 

R3 = 2.*R1*R2 
SI = 1 


C I THIS DO LOOP CALCULATES THE THETA(N) | 

C IAND THE IMAGE SLANT RANGES R8(N) AND R9(N)| 


DO 30 N = 1, N1 

IF(Sl.GT.O) T1(N) = (N-1)*B+G 
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IF(Sl.LT.O) T1(N) = N*B-G 
SI = - SI 

R8(N) = SQRT(D2-R3*COS(T1(N)-D)) 
R9(N) = SQRT(D2-R3*C0S(T1<N)+D)) 
30 CONTINUE 


C I SUM THE PRESSURE OVER ALL IMAGES| 


PI = 0.0 
P2 = 0.0 

DO 40 N = 1, N1 

S2 = (-1)**(INT(N/2)) 

C I-------_________-__| 

C I REFLECTION COEFFICIENTS ALONG NTH UPPER PATH| 


W1 = 2*C2*AL 
II = INT((N-l)/2) 

DO 50 I = 1,11 

S(I) = (R1*SIN(T1 (N)-2*I*B) 

& ^ +R2*SIN(2*I*B-D))/R8(N) 

IF (S (I) .GE.l) S (I)=l 
C(I) = SQRT(1-(S(I)*S(I))) 

T = S(I)/D1 

WO = (-C2+(C(I) *C(I) ) ) 

Y = DSQRT( (W0*W0) + (W1*W1) ) 

Z = WO 
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IF(Y.LE.Z) Y = Z 
Y1 = Q1*SQRT(Y+WO) 

Y2 = -Q1*SQRT(Y-WO) 

Z1 = T-Y2 
Z2 = -Y1 

Z3 = Zl/(Z1*Z1+Z2*Z2) 
Z4 = -Z2/(Z1*Z1+Z2*Z2) 
Zl = T+Y2 
Z2 = Y1 

Z5 = Z1*Z3-Z2*Z4 

Z6 = Z1*Z4+Z2*Z3 
E(I) = Z5 
F(I) = Z6 


50 CONTINUE 

C I —------------- I 

C PRODUCT OF REFLECTION COEFFICIENTS ALONG NTH 

C UPPER PATH 


Zl = 0 
Z2 = 0 
Z3 = 0 
Z4 = 0 
Z5 = 1 
Z6 = 0 

IF(N.LE.2.00) GO TO 110 
DO 60 I = 1, II 
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Z1 = E(I) 

Z2 = F(I) 

Z3 = Z5 
Z4 == Z6 

Z5 * Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
60 CONTINUE 

110 Z1 = Z5 

Z2 = Z6 
T = T4*R8(N) 

Z3 = COS(T) 

Z4 = -SIN(T) 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
PI = P1+S2*Z5/R8(N) 
P2 = P2+S2*Z6/R8(N) 
11 = 11+1 


C [REFLECTION COEFFICIENTS ALONG NTH LOWER PATH| 


DO 70 I = 1, II 

S(I) = (R1*SIN(T1(N)-2*(I-1)*B) 

& + R2*SIN(2*(I-1)*B+D))/R9(N) 

IF(S (I) .GT.l) S(I)=1 
C(I) = SQRT(1-S(I)*S(I)) 

T = S(I)/D1 
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wo = -C2+C(I) *C(I) 

Y = DSQRT{(WO*WO)+(W1*W1)) 
Z = WO 


IF (Y 

.LE.Z) Y = Z 

Y1 = 

Q1*SQRT(Y+WO) 

Y2 = 

-Q1*SQRT(Y-WO) 

Z1 = 

T-Y2 

Z2 = 

-Y1 

Z3 = 

Zl/ (Z1*Z1+Z2*Z2) 

Z4 = 

- Z2/(Z1*Z1+Z2*Z2 

Z1 = 

T+Y2 

Z2 = 


Z5 = 

Z1*Z3-Z2*Z4 

Z6 = 

Z1*Z4+Z2*Z3 

E(I) 

= Z5 

F(I) 

= Z6 


70 CONTINUE 

C I —---------------I 

C PRODUCT OF REFLECTION COEFFICIENTS ALONG NTH LOWER 

C PATH 


Z1 = 0 
Z2 = 0 
Z3 = 0 
Z4 = 0 
Z5 = 1 
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Z6 = 0 


DO 80 

I 

= 1,11 

Z1 

= 

E(I) 

Z2 

,= 

F(I) 

Z3 

= 

Z5 

Z4 

= 

Z6 

Z5 


Z1*Z3-Z2*Z4 

Z6 

= 

Z1*Z4+Z2*Z3 

CONTINUE 


Z1 = Z5 
Z2 = Z6 
T = T4*R9(N) 

Z3 = COS(T) 

Z4 = -SIN(T) 

Z5 = Z1*Z3-Z2*Z4 
Z6 = Z1*Z4+Z2*Z3 
PI = P1+S2*Z5/R9(N) 

P2 = P2+S2*Z6/R9(N) 

40 CONTINUE 

PZ = R1*SQRT(P1*P1+P2*P2) 

PH = ATAN2(P2,P1)*T6 
D = (B-D)*T6 

WRITE(6,3) D, PZ, PH 

3 FORMAT(3X,"D = ' , F8.4, 3X, ' PZ = , F9.5,3X, ' PH = '',F9.5) 
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CONTINUE 


END 



PROGRAM DOUBLTIC 


APPENDIX "C" 


Q*********************************************************** 

C THIS PROGRAM PRODUCES THE OUTPUTS OF THE DOUBLET ANALYSIS 

C OF THE IMAGE THEORY. THE SOURCE AND ITS FIRST SURFACE 

C REFLECTED IMAGE ARE TREATED AS A SEPARATE DOUBLET, AND 
C THE REST OF THE IMAGES ARE TREATED AS PAIRS OF ACOUSTIC 
C DOUBLETS. THIS PROGRAM IS INITIALY MADE TO CALCULATE THE 
C REFLECTION COEFFICIENTS FOR A SLOW BOTTOM WITHOUT ANY 
C ABSORPTION IN THE BOTTOM. THAT WAY, THE REFLECTION 
C COEFFICIENTS ARE REAL, AND THE PROGRAM IS SIMPLER TO 

C INITIALY VERIFY THE VALIDITY OF THIS NEW ANALYSIS. IT 

C ALSO WORKS FOR DOWNSLOPE/UPSLOPE CASE ONLY. THE INITIAL 
C SETUP (BEFORE REACHING THE DOLOOP #20) IS APPLICABLE TO 

C ALL OTHER CASES (CROSS SLOPE, FAST BOTTOM, ABSORPTION IN 

C THE BOTTOM), BUT THE PROGRAM WOULD REQUIRE MODIFICATIONS 
C IN DOLOOP #20 AND BEYOND 

Q********i#r****** * * * * * ******* * * ****************************** 

INTEGER A,N1,NQ1,M1,K,M,MQ,N 

DOUBLE PRECISION B,G,D1,Cl,C2,YO,R1,R2,R3,CC,PI,D2, 

*TQQ,TQQl,PQ2(5000),TQQ2,TQQ3,T4,T6,ROP,SNSOP,P1,P2, 

*Z1(90),Z2(90),GNP(90),THP(90),EP(90),PHIP(90), 
*PSIP(90) ,Z3(90) ,Z4(90) ,T11(90) , RL(90) ,RU(90) ,SL,SU, 
*RLD(90) ,RUD(90) ,DRU(90) ,DRL(90) ,SUMRU(90) , SUMRL(90), 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


*Z5(90) ,Z6(90) ,Z7(90) ,Z8(90) ,RN0P(90) , D (5000) , PQl (5000) 
PI = ACOS(-l.ODOO) 

WRITE(*,*)'INPUT B,G,R1,R2, Y0,D1,C1,C2,A' 

READ (5, *)B,G,R1,R2, Y0,D1,C1,C2,A 
WRITE (6,2)B,G,R1,R2, Y0,D1,C1,C2,A 
2 FORMAT(9X,'B =',F5.2,/, 6X,'G = ',F6.2,/, 6X,'R1= 

*',F5.2,/, 6X,'R2 = ',F5.2,/, 6X,'Y0 = ',F3.1,/, 

*6X,'D1 = ',F4.2,/, 6X,'C1 = ',F7.2,/, 6X,'C2 = 

*',F10.4,/, 6X,/,9X,'A = ' , 14) 

*INPUT PARAMETERS 

FR = FREQUENCY 
X = SCALING DISTANCE 
B = WEDGE ANGLE(DEG) 

G = SOURCE ANGLE MEASURED FROM THE SURFACE (DEG) 

D = RECEIVER ANGLE MEASURED FROM THE SURFACE (DEG) 
N1 = NUMBER OF MONOPOLES(IMAGE POINTS) 

NQl = NUMBER OF PAIRS OF ACOUSTIC DOUBLETS 
R1 = SOURCE DISTANCE(IN DUMP DISTANCES) FROM APEX 
R2 = RECEIVER DISTANCE(IN DUMP DISTANCES) FROM APEX 
YO = SOURCE-RECEIVER HORIZONTAL DISTANCE IN DUMP 
DISTANCES ALONG APEX 
D1 = RH01/RH02 

CC = C1/C2 

A = NUMBER OF RECEIVER POSITIONS 
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C CHOOSE SLOW OR FAST BOTTOM BY VALUE OF SPEED RATIO 

C CC = C1/C2 

Q********************************** ************************* 

Q****:*c)«c****^***********i*c***:«r******************^*******^****’*f 

C MAIN PROGRAM 

Q * vie *** ^ i»r ***** :*c * jic/c * o*r ************** i*c ************* >«c ** 5*c * * 

K = 0 

CC = C1/C2 

N1 = INT(180./B) 

NQl = INT((INT(360./B)-2)/4) 

T6 = 180./PI 
B = B/T6 
G = G/T6 
D2 = 1./D1 
TQQ = TAN(B) 

C 

C ******************************************************** 

C DECISION ABOUT SLOW OR FAST BOTTOM 

C T4 = PI/(2*SIN(ACOS(CC))*TAN(B)) FOR FAST BOTTOM 

C T4 = PI/(2*TAN(ACOS(1/CC))*TAN(B)) FOR SLOW BOTTOM 

Q ******************************************************** 

IF(CC.LT.l) THEN 
TQQl = ACOS(CC) 

TQQ2 = SIN(TQQl) 

ELSE 
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TQQl = AC0S(1/CC) 

TQQ2 = TAN(TQQl) 

ENDIF 

TQQ3 = 2.*TQQ2*TQQ 
T4 * PI/TQQ3 
PRINT*,'KIX = ',T4 
C X = T4*C1/(2*PI*FR) 

CALL EXCMSC'FILEDEF 79 CLEAR' ) 

DATA NOA/79/ 

CALL EXCMS('FILEDEF 79 DISK IFD0UT4 PLOTTER Al') 

R3 = 2*R1*R2 
DO 20 Ml = 1,A 
K = K+1 

C START COUNTING THE RECEIVERS POSITIONS 

D(M1) = (B/A)*M1 
DO N = 1, N1 

C CALCULATION OF THE ANGLE OF EACH IMAGE FROM THE BOTTOM 
IF( (N/2) .LT. (N/2.) ) THEN 
Til(N) = N*B-G 
ELSE 

Tll(N) = (N-1)*B+G 
ENDIF 

C CALCULATION OF THE IMAGE SLANT RANGES 
C LOWER IMAGES 

RLD(N) = DSQRT(R1**2+R2**2+Y0**2- 
& -R3*COS (Tll (N)+B-D(M1) ) ) 
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C UPPER IMAGES 


RUD(N) = DSQRT(Rl**2+R2**2+yO**2- 
& -R3*C0S(Til(N)-B+D(Ml))) 

ENDDO 

C CALCULATION OF THE REFLECTION COEFFICIENTS FOR A SLOW 
BOTTOM WITHOUT ANY ABSORPTION 
C LOWER IMAGES 

DO N = 1,N1 
RL (N) = 1. 

MQ = INT((N-1)/2) 

DO M = 0,MQ 

SKI = R2*DSIN ( (2*M+1) *B-D(M1) ) 

SL = DASIN((R1*DSIN(T11(N)-2*M*B)+SK1)/RLD(N)) 

A50 = DSQRT(CC**2-(COS(SL))**2) 

VL == (D2*DSIN(SL)-A50) / (D2*DSIN(SL)+A50) 

RL(N) = RL(N)*VL 
ENDDO 
ENDDO 

C UPPER IMAGES 

DO N = 3,N1 

C THIS DOLOOP STARTS FROM THE THIRD UPPER IMAGE BECAUSE THE 

C FIRST AND THE SECOND UPPER IMAGES ARE CONSIDERED AS A 

C NEUTRAL DOUBLET 
RU(N) = 1. 

MQ = INT((N-1)/2) 

DO M = 1,MQ 
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SK2 = R2*DSIN((2*M-1)*B+D(M1)) 

SU = DASIN((R1*DSIN(T11(N)-2*M*B)+SK2)/RUD(N)) 
B50 = DSQRT(CC**2-(C0S(SU))**2) 

VU = (D2*DSIN(SU)-B50)/(D2*DSIN(SU)+B50) 

RU (N) = RU (N) *VU 
ENDDO 
ENDDO 

SUM AND DIFFERENCE OF THE REFLECTION COEFFICIENTS 
BETWEEN THE TWO IMAGES OF EACH ACOUSTIC DOUBLET 
DO N = 1,NQ1 

IF ((2*(N+2)).LE.Nl) THEN 
DRU(N) = RU(2*N+1)-RU(2*N+2) 

SUMRU(N) = RU(2*N+l)+RU(2*N+2) 

DRL(N) = RL(2*N-1)-RL(2*N) 

SUMRL(N) - RL(2*N-1)+RL(2*N) 

ELSE 

DRU(N) = 0 
SUMRU(N) = 0 
DRL(N) =0 
SUMRL(N) = 0 
ENDIF 
ENDDO 

PRESSURE FIELD FROM THE NEUTRAL DOUBLET (EQ.31) 
ROP = SQRT{R1**2+R2**2+YO**2-2*R1*R2*COS(D(M1))) 
SNSOP = {R2/R0P)*D(M1) 

PI = (2/ROP)*SIN(T4*R1*G*SNS0P)*(SIN(T4*R0P)) 




J 


c 

c 

c 


c 


c 


c 


c 


c 


c 


P2 = (2/ROP)*SIN(T4*R1*G*SNS0P)*(COS(T4*R0P)) 
CALCULATION OF THE TOTAL PRESSURE FIELD. ADD THE 
FIELD OVER ALL DOUBLETS ACCORDING TO EQS.51 AND 52 
DO 30 N = 1,NQ1 
EQ.35 

RNOP(N) = SQRT(Rl**2+R2**2+Y0**2- 
& -2*R1*R2*C0S(2*N*B)) 

EQ.45 

GNP(N) = T4*R1*R2/RN0P(N) 

EQ.4 6 

THP(N) = GNP(N)*G*SIN(2*N*B) 

EQ.48 


& 


EP(N) = GNP(N)*G*D(Ml)*COS(2*N*B) 

EQ.47 

PHIP(N) = GNP(N)*D(Ml)*SIN(2*N*B) 

EQ.49 

KSIP(N) = 0.5*(D(M1)*EP(N)/G- 

-PHIP(N)**2/(T4*RN0P(N))) 


& 


Sc 


APPLY EQ.52 

Z1(N) = SIN(THP(N)+EP (N))* 

*SIN(T4*RN0P(N)+PHIP(N)+KSIP(N)) 
Z1(N) = SUMRU(N) *Z1 (N) 

Z2 (N) = COS (THP (N) +EP (N) ) * 

*COS (T4*RN0P (N)+PHIP (N)+KSIP (N) ) 
Z2(N) = DRU(N) *Z2 (N) 

Z3 (N) = SIN (THP (N)-EP (N) ) * 
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& *SIN(T4*BN0P(N)-PHIP(N)+KSIP(N)) 

Z3(N) = -SUMRL(N)*Z3(N) 

Z4(N) = COS (THP (N)-EP (N) ) * 

& *COS (T4*RN0P (N) -PHIP (N) +KSIP (N) ) 

Z4(N) * -DRL(N)*Z4(N) 

Z5(N) » SIN(THP(N)+EP(N))* 

& *COS(T4*RNOP(N)+PHIP(N)+KSIP(N)) 

Z5(N) = SUMRU(N) *Z5 (N) 

Z6(N) = COS(THP(N)+EP(N) )* 

& *SIN(T4*KN0P(N) +PHIP(N)+KSIP (N) ) 

Z6(N) = -DRU(N) *Z6 (N) 

Z7 (N) = COS (THP (N) -EP (N) ) * 

& *SIN(T4*RN0P(N)-PHIP(N)+KSIP(N) ) 

Z7 (N) = DFIL (N) *Z7 (N) 

Z8 (N) = SIN (THP (N)-EP (N) ) 

& *COS(T4*RNOP(N)-PHIP(N)+KSIP(N)) 

Z8(N) = -SUMRL(N) *Z8 (N) 

PI = (-1)**N*(1./RN0P(N) )* 

& *(Z1(N)+Z2(N)+Z3(N)+Z4(N))+P1 

P2 = (-1)**N*(1./KN0P(N))* 

& *(Z5(N)+Z6(N)+Z7(N)+Z8(N))+P2 

CONTINUE 
AMPLITUDE 

PQl(Ml) = R1*SQRT(P1**2+P2**2) 

PHASE 


PQ2(M1) * ATAN2 (P2,P1) *T6 






c CONVERT THE RECEIVER ANGLE IN DEGREES 

D(M1) « D (Ml) *180./PI 
C WRITE (6, 3)D (M1) ,PQ1 (Ml) ,PQ2 (Ml) 

C 3 FORMAT(9X,F6.3,8X,F9.6,8X,F10.5) 

20 CONTINUE 

C STORE THE RESULTS IN THE FILE IFDOUT4 PLOTTER A1. FROM 

C THAT FILE THEY ARE RECALLED IN URTEXTIC EXECUTION, FOR 

C COMPARISON WITH THE IMAGE THEORY RESULTS. 

DO 21 J = 1,K 

WRITE (NOA, *)D (J) ,PQ1(J) ,PQ2 (J) 

21 CONTINUE 
END 
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